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[~^ ' The connections between Cosmology and Particle Physics through Dark Matter (DM) have 

received special attention in the last few years for sharpening the physics case of collider 
physics at the TcV frontier. There are many extensions of the Standard Model (SM), 

(~| . which include a new, stable, weakly-interacting massive particle, possibly responsible for 

0^' the observed relic DM in the Universe. The LHC will provide first important data to 

address the question whether one of these scenarios is indeed realised in nature. The ILC 

g^ ' measurements of the properties of a DM candidate and of those other particles participating 

,.J^ , to its interactions in the early Universe may allow us to predict its relic density with an 

accuracy comparable to that currently achieved by CMB observations at satellites. With 
these data in hand, the comparison of the results would have striking consequences for our 

^ ' understanding of dark matter. 

in 

o 

(N 



2 Neutralino Dark Matter Density in MSSM and the ILC 



Supersymmetry emerges as the best motivated theory of new physics beyond the SM. It 
solves a number of problems, intrinsic to the SM and, most important to our discussion, the 
r — ' conservation of R-parity introduces a new stable, weakly interacting particle. The WMAP 

f^ . CMB data, and other astrophysical data, already set rather stringent bounds on Super- 

symmetry parameters, if the neutralino is responsible for saturating the amount of DM 
observed. 

/K^ ' The potential of ILC data at 0.5 TeV and 1.0 TeV for determining the DM relic density, 

^ . ri^, in Supersymmetry has been investigated in detail in [1]. This study selected a set of 

benchmark points, the so-called LCC points, representative of various scenarios and deter- 
mined the 51^ probability density function by a scan of the full MSSM parameter space and 
retaining those points compatible with the measurements available at the LHC and ILC 
within their accuracy. 

3 e+e" -^ H^A^ at LCC-4 with Full Simulation 

We consider here a specific Supersymmetric scenarios, in which the DM candidate is the 
lightest neutralino, Xi ^-nd its relic density is controlled by the rate of neutralino annihilation 
through the CP-even heavy Higgs pole XX ~^ ^- The LCC-4 benchmark point [l] is defined 
in the cMSSM, corresponding to the parameters too=380.00 GeV, TOi/2=420 GeV, tan/3=53, 
A=0, Sgn{^)~+1 and A/fop=178 GeV. We use Isasugra 7.69 [2] to compute the particle 
spectrum and we get M^o=419.4 GeV, M^o = 169.1 GeV and Mf^ = \Qb.h GeV. The e+e" -> 
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H^A° -f bbbb process at ^s = 1 TeV ILC has already been studied for LCC4 [3]. That 
study, based on the parametric detector simulation program Simdet 4.0, showed that the 
A° boson mass can be determined to ±0.8 GeV by imposing the natural width Ta or to 
±2.0 GeV by a simultaneous fit to mass and width. These results, when combined with other 
measurements to be peformed at 0.5 TeV, allow us to predict the neutralino contributoion 
to the dark matter density in the Universe, fi^ to a relative accuracy of 18 % in generic 
MSSM scenarios. 

Here, we repeat the same study on Geant-4-based simulation [4] of the detector response 
and reconstruct the physics objects using processors developed in the Marlin framework [5]. 
This study adopts the LDC detector concept, which employs a large continuous gaseous 
tracker Time Projection Chamber surrounded by a highly granular calorimeter and comple- 
mented by a high resolution Vertex Tracker, for which we have chosen the option based on 
CMOS monolithic pixel sensors. The LDC detector is discussed in details elsewhere[5], the 
design is optimised for achieving excellent parton energy measurements, through the particle 
flow algorithm, and precise extrapolation of particle tracks to their production point. Both 
these features are important to the analysis, which aims to suppress backgrounds by exploit- 
ing its signature 4-6 jet final state and requires good detcrminaton of energy and direction of 
hadronic jets to maximise the resolution on di-jet invariant masses. The jet energy resolution 
has been studied using a simulated sample of single b jets in the energy range from 10 GeV 
to 210 GeV over a polar angle, 0, range 0.4 < 6* < 7r/2, we get SE/E = (0.34± 0.02)/\/l© 
(0.015 ± 0.005), which is consistent with the particle flow performance specifications. Jet 
flavour tagging is based on three observables: the probability for all the particle tracks to 
originate at the event primary vertex, the fraction of the jet energy carried by secondary 
particles and the pt-corrected mass of the secondary particles. These are combined to form 
a discriminant variable which peaks at one for b jets and peaks at zero for non-6 jets. 
At the chosen working point, an 
efficiency for b jets of 85 % 
is obtained with sufficient re- 
jection of lighter quarks to ef- 
fectively suppress the remain- 
ing non-6 background. Signal 
events have been generated with 
Pythia 6.205+Isajet 7.69, in- 
cluding bremsstrahlung effects. 
These events have been passed 
through the full LDC simula- 
tion using the Mokka 06-03 pro- 
gram [7] based on Geant-4. The 
Icio [8] collections produced by 
Mokka have been used as input for 
the Marlin reconstruction. At y^ 
= 1 TeV, the effective e+e^ -^ 
^0^0 production cross section, 
accounting for beamstrahlung, is 
1.4 fb and the decay BR{A° -^ bb) 
is 0.87. The main backgrounds are 
Z^Z", W^W~ production and the inclusive 6666 production. Their cross sections are 0.2 pb. 





Figure 1: Tranverse energy and thrust distributions for 
HA, Z^Z^ and W'^W~ . Generator level distributions 
arc plotted as histograms, results of Mokka + Marlin 
simulation and reconstruction are given for the signal 
process as points with error bars. 
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3.2 fb and 5.1 fb respectively. We assume to operate the ILC at 1 TeV for a total integrated 
luminosity of 2 ab^^. Backgrounds can be significantly suppressed using event shape and 
kinematic variables. We require events to fulfill the following criteria: total recorded energy 
in the event Etot > 850 GeV, total transverse energy Et >350 GeV, charged energy in the 
event Echa > 350 GeV, number or reconstructed particles Ntot >50, number of charged 
particles Ncha >25, event thrust <.95 and Y-^a <0.0025, where ^34 is the 3 to 4 jet crossover 
value of the jet clustering variable. The distributions of some of these variables is shown 
in Figure [1]) for backgrounds and signal, for which a comparison of the generator-level and 
reconstructed values is also given. 

After event selection, the di-jet pairing which minimises the di-jet mass difference has 
been chosen. The di-jct mass resolution has been improved by applying a 4-C fit. We 




500 550 



Figure 2: Dijet invariant mass distribution for e^e~ — + H^A^ events selected by the analysis 
cut. Mass constraint fit and jet flavour tagging has been applied. The distribution for fully 
simulated and reconstructed events (points with error bars) is compared to that obtained 
with parametric simulation (histogram). 

have ported the PUFITC algorithm, developed for the DELPHI experiment at LEP, into 
a dedicated Marlin processor. The algorithm adjusts the momenta of the jets given by 
Pp = e°^pM+bpB+cpc where pV is the fitted momentum, pM is the measured momentum,7?B 
and PC are unit vectors transverse to pM and to each other, and a, b and c are the free 
parameters in the fit. The adjusted momenta satisfy a set of constraints while minimizing 
the fit x^, which is given by S^ (a^ — 00)^/0-^ + fof/o'^ + cf /cr^, where uq is the expected energy 
loss parameter, aa is the energy spread parameter and (Tb,(Jc are the transverse momentum 
spread parameters. In this analysis we use the following constraints: Px ~ Py ~ and 
iJ ± |pz| = y^, where the last condition accounts for beamstrahlung along the beam axis, 
z. We report here preliminary results from the analysis of a sample of 1050 fully simulated 
signal events. After applying final selection and mass constrained fit, the sample of events 
in the region 150 GeV< Mjj <550 GeV gives a selection efficiency of 23 % for signal bbb 
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decays. The resulting mass distribution is shown in Figure [2 We describe the signal as 
a CRYSTAL BALL (CB) function [5] and extract the A° mass, Ma, and width, Ta have 
been by a multi-parameter fit leaving the CB parameters free. We determine the A° mass 
as (419.1±0.9) GeV. This result is remarkably close to that obtained in the earlier analysis, 
based on parametric detector simulation. The production and analysis of fully simulated 
and reconstructed background samples is currently under way. 
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4 Further Constraints on Q^ 

The constraints on LCC4 derived from mass measurements at the LHC and ILC, provide a 
prediction of the DM density in the Universe to a relative accuracy of 18 % with a generic 
MSSM 

This accuracy is still far from that 
achieved by CMB study with satel- 
lites. The main contribution to 
the remaining uncertainty is the 
weak constrain which the data pro- 
vide to MSSM solutions where 
r2j(. is significantly lower than its 
reference for LCC4. A detailed 
study shows that these solutions 
are all characterised by large val- 
ues of the stau trilinear coupling, 
Atau- In the MSSM the f cou- 
pling to the i7° and A'^ bosons 
scales as Atau ""l", + M^ and 

•^^^ cos beta ^cos/3 

Ata„ tan /3 -I- ^, respectively. In the 
funnel region the main annihila- 
tion mechanism is x'^x" —^A'^^ 
bb and Ma < Mf^ +Mf^ . The only 
A^ decay into fs, allowed by CP 
symmetry, A^ -^ fif2 is kinemati- 
cally forbidden. At large values of 
l^taiil, the stau decay process through the H^ 

This channel contributes to the neutralino annihilation rate through x^x^ -^ H'^' -^ Tifi, 
bringing down the corresponding relic density, as observed in the MSSM scans. At the same 
time, a dctermnation of the branching fraction of the decay H'~' —> fifi, allows to constrain 
|Aj(jij|. Figure [3] shows the decay branching fractions of the A" and iJ" bosons computed us- 
ing the HDECAY 2 . program [TU] as a function of Atau- Now, a large H'^ — *■ fifi — > tx^tx'^ 
yield can be detected by a standard bbrr analysis. A preliminary study shows that the A^, 
H^ -^ TT branching fraction can be determined to ± 15 % and A", H'^ -^ bb to zL 7 %, 
from which a limit l^taul < 250 GeV can be derived. This constrain removes the tail at 
low values of O^ and results in a prediction of the neutralino relic density with a relative 
accuracy of 8 %. A detailed study on full simulation to support these preliminary results is 
currently under way. 
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Figure 3: i/° and A^ decay branching fractions as a 
function of the stau trilinear coupling Atau predicted 
by HDECAY. AU the other MSSM parameters have been 
kept fixed to those corresponding to the LCC4 point. 



fifi gets a sizeable branching fraction. 
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